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I. INTRODUCTION 

One of the most successful models for the origin of the primordial energy density 
fluctuations which gave rise to the observed structures in the universe is given by 
inflation. Quantum fluctuations of scalar fields present during the inflationary era 
give rise to classical fluctuations in the energy density. In the simplest model of 
one inflaton field driving the universe expansion, as it slowly rolls down a smooth 
potential, the resulting fluctuations are adiabatic, with scale invariant spectrum and 
Gaussian distributed. These are also the kind of primordial fluctuations usually 
considered in one of the most accepted scenarios for galaxy formation, the standard 
Cold Dark Matter (hereafter CDM) one. This scenario successfully explains many of 
the observed properties of galaxy clustering. However, it has shown to have difficulties 
in explaining some large-scale observations such as the galaxy peculiar motions , 1 the 
clustering of clusters , 3 the angular correlation function of galaxies 3 and some peculiar 

features in the very large-scale galaxy distribution . 4 

Different models in which the inflationary paradigm can be realized have been pro- 
posed. These are generally identified as old inflation, new inflation, chaotic inflation 
and extended inflation. In the old inflation picture the universe undergoes a first order 
phase transition during which the regions in the false vacuum phase expand exponen- 
tially, it has the problem that, due to the speed of the false vacuum expansion, the 
true vacuum bubbles fail to percolate. In the new and chaotic inflationary scenarios 
this problem is solved, as inflation occurs during the slow rolling of the inflaton field 
to the minimum of its potential. The extended inflation model has recently been pro- 
posed as an alternative solution to the problems of old inflation . 5 A first order phase 
transition for the inflaton field is considered in the frame of the Brans-Dicke theory of 
gravity. This makes the universe expand as a power law (instead than exponentially) 
during the phase transition, so that the percolation of false vacuum bubbles can oc- 
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cur. However, this model has the problem that it is not possible to simultaneously 
satisfy the bounds coming from the observed smoothness of the Cosmic Microwave 
Background Radiation (hereafter CMBR) and from time-delay experiments. Modified 
versions of this model have been proposed to overcome this difficulty. 8 

Later on it has been pointed out by Linde 7 that the Brans-Dicke gravity theory also 
has interesting consequences for chaotic inflationary scenarios due to the fact that 
the effective gravitational constant can take different values in different places of the 
universe. One of the modifications of the original extended inflation model proposed® 
considers more general scalar-tensor theories of gravity in which the Einstein gravity 
theory is an attractor for the solutions, but that can considerably differ from it in 
the past. The main purpose of this paper is to explore the consequences of chaotic 
type inflation taking place in these gravitational theories. We show that the fact that 
curvature fluctuations are determined by both the fluctuations in the inflaton and in 
the Brans-Dicke field gives rise to perturbations with non scale-invariant spectrum. 
In particular, it is possible that the spectrum has a peak at a given wavelength. This 
kind of spectrum can help to reconcile CDM predictions with large-scale observations 
without running into conflict with the CMBR anisotropy limits, provided that the 
peak in the spectrum appears at an appropriate wavelength. Other models in which 
this kind of spectrum can arise have been discussed in the literature. 9 For example, a 
mountain on the top of an underlying scale invariant spectrum arises when there is a 
hill on the scalar field potential. 10 Another possibility is to consider two interacting 
scalar fields. In this case the interaction between the fields can give rise to variations 
in the effective masses of the fields during their evolution and can make change 
sign, thus originating a peak in the resulting density perturbations. (Other kinds of 
non scale-invariant fluctuations in inflationary models have been discussed in Refs. 
11.) The mechanism considered here is quite different. The peak in the spectrum 
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originates when the effect of the inflaton fluctuations becomes dominant over the 
Brans-Dicke field fluctuations in determining the total curvature fluctuations. A 
problem of these models is that, in order to make the features in the spectrum appear 
at scales of astrophysical interest, it is necessary to impose rather precise conditions on 
the inflationary model, either on the initial conditions of the fields or on the potential 
parameters. This problem is also present in the model discussed here, as the initial 
value of the Brans-Dicke field determines the scale at which the peak in the spectrum 
appears. Another model where non scale-invariant perturbation are generated has 
been recently proposed, 13 in which a period of extended inflation is followed by a 

period of slowroll inflation. 

The organization of the paper is as follows. In section 2 the chaotic extended in- 
flation model is introduced; the conformal transformation to the Einstein frame, the 
equations of motion in that frame and the inflationary solutions are discussed. In 
section 3 the generation of curvature perturbations in this model is studied and the 
possibility that their spectrum is non scale-invariant at observable scales is discussed. 
The resulting spectrum of perturbations for some particular values of the model pa- 
rameters is presented. In section 4 the evolution of the inflaton and of the Brans-Dicke 
field is analysed in the frame of the stochastic approach to inflation. In section 5 we 

discuss the results. 

II. CHAOTIC EXTENDED INFLATION MODEL 
A. The model 

We wiU explore the realization of chaotic inflation in the tame of a class of scalar- 
tensor theories of gravity which are described by the action 

S = J + \d*todnr - V(<r)^ , (1) 
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where <f> is the Brans-Dicke field and <r the inflaton field. The scalar field <f> is nor- 
malized to be (2/C 2 )- 1 at present, with K 2 = 8wG N (Gs being the Newton constant). 
In the standard Brans-Dicke theory w($) is a constant, and it is constrained by 
present observations to be large, u > 500. A more general scalar-tensor theory is the 
Bergmann 13 -Wagoner 14 -Nordtvedt 15 one, where u((f>) is an arbitrary function. Special 
cases are Bekenstein’s 16 “variable mass theory” and Barker’s 17 “constant Gft theory”. 
In these theories the function w(<£) can have the property that, for the present value 
of the scalar field <j> } the theory is indistinguishable from General Relativity, but, 
for past values of it could lead to significant differences in cosmological models. 
Particularly interesting is the case in which the system dynamically evolves to an 
attractor solution indistinguishable from General Relativity. This kind of theories 
has recently been considered in order to overcome the problems associated with the 
original extended inflationary model proposed by La and Steinhardt. 8 With this scope 
a particular class of models, based on the function w($) defined by 

+ 2 = (1 -2/C^)**’ ^ 

with /? and fi constants, has been analysed by Garaa-Bellido and Quiros. 8 We will 
analyse the consequences of chaotic inflation in the frame of this gravity theory. 

The action in eq. (1) is written in the Jordan frame, where the gravitational 
coupling depends on <f>. However, it shows to be more convenient to work in the 
conformally transformed Einstein frame in which the gravitational part of the action 
takes the usual form. The conformal transformation connecting the two frames is 
given by 

g a b = fl -a (ja4, with D = V^- ( 3 ) 
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It is also useful to redefine the scalar field as 18 


In this frame the action (1) transforms to 

s = / + \r h s.f8, f + - sr»vw) . (5) 

We will call F{ip) = = (2£ 3 $) *• 


B. Evolution of the fields 

The equations of motion for the Brans-Dicke field <p and the inflaton cr in the 

Einstein frame are given by 

<r 3 

<fi + ZH<p = -2F(<p)d v F V(ff) + a/y, (6) 

* + 3 H* + j*=-F{'p)d.V, (7) 

having restricted our analysis to a homogeneous and isotropic umverse. Here a dot 
denotes differentiation with respect to t and H = a/a. The Friedmann equation reads 

H 3 = y + F(v>)y + • ( 8 ) 

To deal with this coupled system of equations it is convenient to change the time 
variable from t to a = ln(a/a 0 ), where a 0 can be taken as the value of the scale factor 
a at the beginning of inflation. The system transforms to 

c T t 2 

H 3 <p" + K 3 F 3 V<p' = -2 Fd v F V + H 3 d v F — , (9) 

ffV + K’F’W + JI'MyV = -F(v)d,V, (10) 

and 

_ (ii) 

_ 3-f (*>» + F(*>K)’ 
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where primes denote derivatives with respect to a. 

In what follows we will consider the particular function w(<f>) quoted in equation (2). 
Two particular constant fi values make it possible to handle the problem analytically: 
these are /i = 1 and \l = 2. The corresponding F(<p) is given by 

i r M=i(v ? ) = coshS ( x: ¥ J V^)’ ( 12 ) 

F^M = 1 + e-Ww. (13) 

In the first case the field <p rolls down its potential and asymptotically approaches 
ip = 0 at late times. In the second case the field rolls down asymptotically approaching 
infinity. All the results obtained are qualitatively similar in the two models, so for 
definiteness we will develop the computations for the p, = 2 case. 

We will consider for the inflaton a potential V((r) = Af 4 exp(— 7/C<r), with 7 a 
positive constant, which gives interesting possibilities for the spectrum of the resulting 
energy density perturbations. 19 

It can be seen that the condition to have an inflationary expansion (o > 0) corre- 
sponds to 

2 > K\<p n + F &)**). (14) 

In this regime we can neglect the second term in the right hand side of eq. (9) and 
the second derivative term, so that the equation decouples from the one for <r and 
can be integrated, leading to 

~ <Pi) + exp (ipKy/yp) - exp(<p f Ky/2/0) = |(a - a,). (15) 

It can be seen that neglecting the second derivative term is a good approximation 
for 4/3 <C /?( 1 + exp(Ktp ^2//?)) J and that neglecting the inflaton kinetic term is a 
good approximation for 7 3 4(1 + e ~ Klfi Note that for large negative values 
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of p eq. (15) reduces to the solution corresponding to a scalar field rolling down an 

exponential potential, i.e. p = pj + 2-j2/(3K 2 (a - <*/)• 

At late times p > 1 and F(p) approaches unity. The effective Newton constant 

tends to its present value as 0 - (2 K')' 1 and from eq. (2) we see that u, grows up 
making the theory indistinguishable from General Relativity. 

The evolution equation for a can be solved in the slow-rolling approximation ne- 
glecting the second derivative term. It is convenient to solve it in terms of p instead 
of a. Combining eqs. (10) and (9), in the approximation discussed above, we get 


dp d^F (i6) 

da 7JC ’ 


which can be directly integrated. The solution is given by 

a = a, + ?£ (e K *^ - • ( 17 ) 

The approximations performed hold for 8/3 < 0(1 + exp(/C^^2/0)) and 7 < 

6(1 -I- cxp(—Ktp^2/ft))' 

As we will see in the following, the relevant quantities for the study of density 
perturbations are p n and Fa n , which in this case are given by 




£ 3 0( l + e^V 5 ^) 3 ’ 


Fa n - 


/C 3 (l + e~ Kv ^) 


(18) 

(19) 


From eqs. (18) and (19) we see that p n is a decreasing function of p (and hence of 
time), while Fa 17 is an increasing one. Thus, for early times p n will be the dominant 

one, while Fa 12 will dominate at late times. 

The exponential potential for the inflaton, as opposed to other chaotic inflation 
potentials, does not lead to a phase where the inflaton field oscillates starting the 
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reheating process. The end of inflation in this case must be determined by some ex- 
ternal reheating mechanism. Without loss of generality we can consider that inflation 
ends at 07 = 0, in which case M 4 corresponds to the potential energy at the end 
of inflation. On the other hand, we know that at the end of inflation the value of 
the Brans-Dicke field <f> in our patch of the universe must be such that the gravity 
theory is now very close to General Relativity. As the strongest variation of the <p 
field occurs during the inflationary era, we will impose this constraint at the end of 
inflation, F(<pf) « 1. 

To conclude this section let us briefly discuss the relation between quantities in the 
Einstein frame and in the Jordan (physical) frame. If we write the line element in 
the Einstein frame as 

da 1 = dt 3 - a\t)dl 3 
and the one in the Jordan frame as 

da 3 = dr 3 - b 3 (r)dt 3 , 

they are related by 


= r). 

(20) 

il = t. 

(21) 


At late times, when F approaches unity, the two frames practically coincide. 

From eq. (20), we see that the number of e-foldings of inflation in the two frames 
is related by 

In .(«/«o) = \ WFo/F f ) + ln(6//6 0 ), (22) 

where quantities at the beginning of inflation have been denoted by a subscript 0 and 
quantities at the end of inflation by an /. 
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Furthermore, the condition for inflation in the Einstein frame (d’u/dl 1 > 0) is 
fulfilled when condition (14) holds. It can be seen that under this condition also the 

Jordan frame inflates (d 5 6/dr l > 0). 


C. Constraints from Post-Newtonian experiments 
There are different constraints that an alternative theory to General Relativity must 
fulfill in order to agree with post-Newtonian experiments. 20 They constrain essentially 
the present value of u; and the present rate of change of the Newton constant. 

As the theory that we are considering evolves approaching General Relativity, if we 
impose these constraints at the end of the inflationary era, they will be also fulfilled 
at the present time. In the last subsection we have seen that the theory approaches 
General Relativity during inflation. It is easy to see that this is also the case during 
the radiation and matter dominated eras. The equation of motion for <p during the 
radiation dominated period can be written in the Einstein frame as 

1 (aV) = -^(/>-3p) = °, ( 23 ) 

where we have denoted d/dt by a dot. Changing the time variable to a, this equation 
can be solved as 

y , = y >, + /,(1 - ( 24 ) 


where <p f and denote the Brans-Dicke field and its derivative at the end of inflation. 
As v' t is positive, we see that y> continues to increase during the radiation dominated 
era and consequently F(<p) - 1, * - (2K 3 )' 1 and *-<*>. It can be seen that p 
also grows during the matter dominated era. The equation of motion (23) in this case 

can be written as 


3 , h/VP 

+ 2 V 2 K i + 


(25) 
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which is difficult to integrate; from it, however, it is easy to see that <p' is always 
positive as 


ip'e 3 *' 3 = <p' r e 3ar/3 + 


2K 


[“da'— 

J* r i 


»3a'/J 




(26) 


where the subscript r denotes the values at the end of the radiation dominated era. 
The last term is always positive, thus ip' > <p' r e 3 ( aT ~ a V 3 > 0. (For a demonstration of 
the convergence to General Relativity starting from the Jordan frame, see Ref. 8.) 

The first constraint to be considered is the one coming from the present value of u 
determined from radar time-delay measurements, u > 500. 31 The value of u> at the 
end of inflation is determined through the relation 


/u> + 3/2 

~~T~ 


= 1 + e 




Thus, the constraint on u> determines a lower limit on the final value of the Brans- 
Dicke field at the end of inflation, corresponding to 

> In - l) . (27) 

Another constraint comes from the present rate of change of the Newton constant. 
The observable gravitational constant is given in terms of <f> by 20 




2uj 3 


)• 


(28) 


Its rate of change per expansion time in the Jordan (observable) frame is given by 


dG 1 yj' d v F f 2/3 + 1- 1/F 3 \ _ 
dr G H ~ F3/a ^20 + (1-1 /F)*) =V 


(29) 


Experimental limits coming from the Viking lander ranging data 32 constrain it to 
ij = 0.04 + 0.08 (assuming H 0 = (2 x 10 lo yr) -1 )- From the analysis of the <p evolution, 
we see that tj evolves towards zero in all the periods of interest. Thus, if at the end of 
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inflation r, was inside that interval, it will be there even more now. In order to fulfill 
this condition at the end of inflation it is necessary that 

Kv,fii0 > 1" - l) . < 30 > 

which «lw.ys holds if condition (27) is satisfied. Big Bang nucleosynthesis imposes 
a stronger limit on the variability of (5, 11 |,| < ±0.018. This bonnd amounts to 
changing the factor 10 by a factor 15 in eq. (30); this condition is also fulfilled 

whenever condition (27) is. 

Finally, let us consider the constraint coming from the lower limit on where 

H 0 is the present Hubble constant and r 0 the cosmic time, Hfo > 0-4- Making use 
of the Einstein equations 



it can be seen that, for values of <p at the end of inflation like those allowed by eq. 
(27), this constraint is easily fulfilled. 

In the next section the generation of density perturbations is studied working in the 
Einstein frame, in which the treatment is closely related to the usual chaotic inflation 
case.* 4 The fact that the two frames coincide at late times, when all the wavelengths 
of interest enter the Hubble radius, implies that our results apply without changes to 
the observable perturbations in the Jordan frame. An alternative method to study 
density perturbations has been developed, in which the Jordan frame is used.” 

III. generation of perturbations 

A. Curvature perturbation* 
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In this model the energy density perturbations arise as a consequence of fluctuations 
in the inflaton field and in the Brans-Dicke field. 

We will work in the longitudinal gauge where the metric takes the form 

da 7 = (1 - 2$)dt 7 + a J (t)(l + 2$)dl 7 , (32) 


where $ coincides with the gauge invariant potential defined by Bardeen 58 and 
corresponds to the usual gravitational potential. It can be expanded in a superposition 
of modes of wavenumber k, $£ (we will not include the subindex k in what follows). 

Curvature perturbations are described by the amplitude of the variable £, 5r given 
by 

Ht7^ 4+ *''*>W i+ 5(T^W’ < 33 > 

which is constant during the evolution for wavelengths much larger than the Hub- 
ble radius; w is the ratio between the pressure and the energy density. The power 
spectrum of £ (variance per In k ) is defined as 

P<W = ICkl* > • (34) 


For scales that are just entering the Hubble radius during the matter dominated 
era the amplitude of P( and P# are related by (see, e.g., Refs. 27,10) 


P 1 / 2 

U 


HC 


5 

3 


” pi / 2 

— - 


HC 


5 6p 

2 7 


HC 


15 ST 

~2 T 


(35) 


HC 


where the subscript HC refers to the Hubble radius crossing time. The amplitude of 
£ modes during inflation in terms of the perturbations of <p and a can be computed 
as follows. From the 0 — j component of the Einstein equations we have 


• K 7 

$ + H$ = — — {<p8<p + F{tp)crS<r) . (36) 
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Thu., a good estimate of < for large wavelengths during the inflationary period is 
given by 


? i , • . H(<p6<p + F(<p)<r6<r) 

<“!rh; ( * +jr #) “ — 


(37) 


or in terms of the variable a 

(/*¥> + FjipySv) (38) 

» ” + F ( v ’) <r ' 3 

This treatment is a good approximation in our one, as both the fields and their 
derivatives are smoothly varying during the entire period of inflation. A more careful 
treatment, integrating the evolution agnation for the field, art metric perturbations, 
is instead necessary in other cases. 


B. Power spectrum 

The amplitude of the held fluctuation. t V and «<r is given by the quantum fluctu- 
ation. of the short wavelength modes of these fields. The power spectrum per Ink 
of scalar field fluctuation, can be computed from the two-point function of the Held. 

For a power-law inflation it is given by 

. . i-3v 

, ( 39 ) 


PM = 


k s p irMI 1 

k p/2 

(2ir) 3 (p-l) Ha* 

aH p — 1 


where p is the power of the expansion (a <x 9 and p = 2/y’ for the exponential 
potential V oc exp <-,*,) case) and pm (3p - l)/(2p - 2). For p - oo it reduces to 
the well-known value HV(2*)». corresponding to a d« Sitter space. At Hubble radius 
crossing time the power spectrum is given by 


p m _ ffljO!! ( -F/L-V'”' ff» 

(2w) 3 V(p - 1)J 


(40) 


I HC 


where the Huhble constant is evaluated at the value that V get. when the wavenumber 
k crosses the Hubble radius. This amplitude appUes only to minimally coupled scalar 
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fields. Thus, as it has been argued in Ref. 24, it can be applied to the Brans-Dicke field 
fluctuations 8<p in the Einstein frame but not in the Jordan one. For the fluctuations 
of the inflaton field 6<r we cannot take that amplitude, due to the fact that its kinetic 
term does not have the canonical form, because of the factor F(<p) in its Lagrangian. 
However, for the inflaton field it is easier to work in the Jordan frame, where it is 
an ordinary minimally coupled field. Thus, in this frame we can estimate the power 
spectrum of 8j<t by eq. (39). Transforming to the Einstein frame we get from eq. 
(40) that, at Hubble radius crossing, 



(see also Ref. 29). 30 In the model that we are studying the universe expansion is not 
given by a single power law for the whole evolution but it can be well approximated 
by different power laws for different periods with slightly different values of p. The 
value of p can be determined from the Friedmann equations 1/p = 3(1 + w)/2 = 
K*{tpP + F<r n )/ 2. 

From eqs. (38) and (41) it can be seen that the curvature perturbations will be 
dominated by the inflaton fluctuations in the case that ip 1 < y/Fcr ' and by Brans- 
Dicke field fluctuations in the opposite case. In the last section we have seen that the 
kinetic contribution of the Brans-Dicke field is the largest one at early times while the 
inflaton dominates at later times. Hence, curvature perturbation at larger scales are 
determined by Brans-Dicke field fluctuations ss —Stp/p', while shorter scale ones 
are determined by inflaton field fluctuations « —Sa/a*. 

1. Qualitative analysis 

The qualitative shape of the spectrum at Hubble radius crossing can be determined 
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as follows. For the larger wavelengths 

(1 + c ^v^)exp - e^^)) » 


(42) 


where C(p) = (2 - 2/p)*'- 1 |r(»')| a /* “ * factor that 8 oes to onc for ? > 1 and is 
smaller than one for values of p closer to one. The spectrum (42) is a decreasing 

function of it has a minimum near V equal to zero, then it increases again and 
it eventually becomes decreasing for large values of *>. (For some values of 7 and 0 
another minimum in P < can appear at a negative value of *, however this is not the 
case in the examples that we will consider). The Hubble radius crossing condition for 
increasing wavenumbers k corresponds to decreasing values of*>. Thus, the spectrum 
of perturbations first decreases with k, it has a minimum and then increases again. 
However, when the kinetic term of the inflaton becomes dominant the curvature 


perturbations are no longer determined by eq. (42), but by 

a + e ‘ W57ii ) 


3/3 


(43) 


which is a decreasing function of tp. Thus, for the shortest wavelengths leaving the 
Hubble radius during the last part of inflation, the spectrum of perturbations again 


decreases with k. Combining the two behaviours we see that a peak in the spectrum 
appears at the wavelengths that are leaving the Hubble radius when the inflaton 
kinetic term becomes dominant. The peak will appear provided that the inflaton 
kinetic energy becomes dominant for a value of *> larger than that corresponding to 
the nimmvm of P<, (eq. (42)). In order to make this happen there is a condition to be 
satisfied by the parameters appearing in F(*>) and V'(o), namely fit 1 £ 4/3. On the 
other hand, the condition that the peak in the spectrum appears at astrophysicaUy 
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interesting scales requires a particular initial condition for the Brans-Dicke field, or 
equivalently, a particular final value of (p at the end of inflation. 

For large wavenumber perturbations, which cross the Hubble radius when the infla- 
ton kinetic term dominates, we can explicitly compute the spectrum of £ perturbations 
at Hubble radius crossing as a function of k from eq. (43) and the crossing condition 
k = aH. It is given by 


> 1 / 3 , 




= ( h/1 \ 

2tt\/37 \KM 7 afj 
with oy the value of the scale factor at the end of inflation. 


(44) 


The spectral index of £ is related to the spectral index of the energy density per- 
turbations n by P^ 2 oc fc( n-1 )/ 3 (for n = 1 the spectrum is scale invariant and P^ 2 
is independent of k). Thus, we see that the spectral index corresponding to the long 
wavenumber region is approximately given by n ~ (1 - 37 J /2)/(l - 7 3 /2), which is 
the well-known result for power-law inflation with an exponential potential. 19 

It has been shown in Ref. 31 that a power spectrum with a primordial spectral index 
n smaller than one, like the one we obtain for the large wavenumber perturbations, 
can help to avoid some difficulties of the CDM scenario, as it has more power on 
large scales than the usual scale invariant spectrum. However, lowering n increases 
the amplitude of fluctuations at large scales, running into possible conflict with the 
CMBR anisotropy limits. This problem may be avoided in the model discussed here 
if the peak in the power spectrum occurs at a scale somewhat larger than 100 Mpc 
so that the amplitude of perturbations decreases for the larger scales relevant for 
the quadrupole anisotropy. The amplitude of the perturbations is determined by the 
value of M. It must be fixed with the help of eq. (35) in such a way that the limits 
on the CMBR anisotropy are satisfied. 
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2. Examples 


The shape of the spectrum for the large wavelengths does not have a simple an- 
alytical expression. Fig. 1 shows the power spectrum of C at horizon crossing for 
the case /3 = 4 and q = 0.38. As these values of the parameters are very close to 
the limits under which the analytical solutions of the equations of motion for f and 
u presented in section II-B are no more valid, the exact numerical solution of the 
system of eqs. (9) and (10) has been used. The value of the Brans-Dicke Reid at 
the beginning of the period of inSation relevant for our patch of the universe needs 
to be in the range between -3.5 /K and -1.5 /K in order to make the peak in the 
spectrum appear at interesting scales. The corresponding value of f, is -f,K ~ 5.7, 
which satisfies eq. (27). We have taken M = 3.4 x 10->m, that make, the amplitude 
of perturbations consistent with CMBR anisotropy limits. The value of ^ determines 
the spectral index for the shorter wavelength region of the spectrum. The considered 
value ^ = 0.38 corresponds approximately to n = 0.85. To decrease the amplitude 
of perturbations at larger wavelengths, small values of 0 are preferable. However, in 
order to make the universe have an inflationary expansion for large negative values 
of y>, there is a constraint on 0, i.e. 0 S. 4. When this constraint is fulfilled, as in the 
case considered in Fig. 1, the inflationary expansion can start at the Planck energy. 

On the other hand, when values of /? smaller than 4 are considered, this does not 
mean that inflation does not happen at all, but just that it starts when the field V 
becomes larger than a certain value, making the number of e-foldings smaller. In 
this case, the evolution of f for large and negative values is too fast, making the 
kinetic energy contribution the dominant one and the expansion non-inflationary. As 
y, evolves, its kinetic term decreases and, when it becomes subdominant, inflation 
starts. Fig. 2 shows the spectrum of perturbation, arising for a case like this. The 
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parameters are chosen to be /? = 3, 7 = 0.43. The value of the Brans-Dicke field at 
the beginning of the last 63 e-foldings of inflation needs to be in the range —1.5 /K 
and — 0 . 5 //C. The corresponding value at the end of inflation is given by Ktp/ ~ 5.4, 
which satisfies the bound coming from eq. (27). The value M = 2 x 10 -3 mp satisfies 
CMBR anisotropy limits. The spectral index corresponding to the short wavelength 
part of the spectrum in this case is given by n as 0.8. The smallest wavenumber 
appearing in the graph corresponds to the first wavelength leaving the Babble radius 
when inflation starts. A necessary condition for the model to be consistent with 
observations is that such a wavelength is larger than our present horizon. 


IV. STOCHASTIC EVOLUTION 

The dynamical evolution of the inflaton and of the Brans-Dicke field, taking into 
account the effects of quantum fluctuations is most clearly studied in the frame of the 
stochastic approach to inflation . 33 This is also the best tool to obtain quantitative 
information about the statistics of curvature perturbations . 33 ’ 34 
In this frame the dynamics of the system is described by two coupled Langevin 
equations for long wavelength modes (k << aH ) of <p and cr. These equations can be 
easily written in the Einstein frame by adding to the classical equations of motion a 
noise term whose amplitude is fixed by the r.m.s. fluctuation of the fields at Hubble 
radius crossing, given by eqs. (40) and (41). They read 


<fi' = 


1 jl 1 £M 2 T( 

1C V 0 (i _|_ 27T 3 / 1 


<j = 
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where rfrp and tj# are Gaussian noises with zero mean and correlation function 
( Vv( a )v<p ( a> )) = (v*( a )v = S(a — a'). In the approximation leading to these 
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stochastic equations the fine-grained components of the fields are treated as free even 

when <p and <r interact, thus we assume (i? v (a) vA a ')) = °- 

By integrating these Langevin equations it is possible to compute the statistical dis- 
tribution of the field fluctuations and of the density perturbations. 34 Non-Gaussian 
fluctuations in the energy density are expected to arise whenever the effect of short 
wavelength quantum fluctuations of a field becomes important compared to the clas- 
sical force in the evolution of the long wavelength part of the field. In the case of 
a single scalar field this criterion implies that non-Gaussian fluctuations are always 
associated to large amplitude fluctuations. 38 This can be easily seen by the analog of 
eq. (38) for one field, C * -KHipfip'. In order to have C small enough to be consistent 
with CMBR anisotropy limits (C?(10- 4 )), it is necessary that the change in ip due to 
quantum fluctuations in one expansion time Sip ~ tf/2ir is much smaller than the 
change ip ' due to the classical force in the same interval. This implies that fluctuations 
of reasonable amplitude are Gaussian distributed. A similar argument applies also 
to our case. As it has been discussed above, curvature perturbations are dominated 
by the fluctuations in that field whose kinetic term gives the largest contribution at 
the time when the associated wavelength leaves the Hubble radius. (For the larger 
scales C * -Sipl<p’ and for the shorter ones C « -Sv/a 1 .) For each of these regimes 
the same analysis applies. At the transition between the two regimes (when the peak 
in the spectrum appears) we do not expect non-Gaussian features to appear in the 
curvature perturbation distribution since all variables change smoothly there. These 
qualitative arguments are fully confirmed by the numerical integration of eqs. (45). 
Using the technique described in Ref. 34 we have numerically integrated the system 
of equations for the two choices of parameters corresponding to Fig. 1 and Fig. 2. 
Initial conditions have been chosen such that only fluctuations that are now inside a 
patch of the universe of the observable size are considered. From the distribution of 
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the inflaton and the Brans-Dicke field the resulting distribution of ( was computed 
using eq. (38). In both cases the distribution of curvature perturbations is pretty 
Gaussian. 


V. CONCLUSIONS 

Chaotic inflation in the frame of scalar-tensor theories of gravity has been anal- 
ysed. We have considered a class of models in which the gravity theory has General 
Relativity as an attractor but that differs from it significantly during inflation. In 
particular, the generation of curvature perturbations in these models was studied. 
These perturbations get contributions from both fluctuations of the inflaton and of 
the Brans-Dicke fields. This makes it possible for the power spectrum to be quite dif- 
ferent from the scale-invariant one. The quantum fluctuations of the field that has the 
largest kinetic contribution when a given wavelength leaves the Hubble radius during 
inflation are the dominant ones in determining the amplitude of curvature pertur- 
bations. As a result, the long wavelength part of the spectrum is dominated by the 
Brans-Dicke field fluctuations, while the short wavelength part by the inflaton ones. 
The spectrum of the curvature perturbations associated to the Brans-Dicke field has 
a dip for the models considered. The spectrum associated to the inflaton fluctuations 
instead is always decreasing with k. Thus, it is possible that a combination of these 
two behaviours results in a curvature spectrum that is peaked at a certain wavelength, 
which might help to reconcile the CDM model with observations. In fact, a peak in the 
spectrum is easily obtained. However, the scale at which the peak appears depends 
quite strongly on the value of the Brans-Dicke field at the end of inflation and the 
height and width are determined by the potential parameters. In order to make the 
peak occur at astrophysically interesting scales particular initial values for <p have to 
be chosen. The final value of the Brans-Dicke field is constrained by post-Newtonian 
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experiments as shown by eq. (27). This constraint is automatically fulfilled by the 
values of <p f coming from the initial values required to have the peak at an interesting 
scale. However, apart from that constraint there is no other restriction on the value 
of w Thus, the reason why the initial value for the Brans-Dicke field should be in 
the allowed range cannot be explained within the model. Finally, we have analysed 
the evolution of the system in the stochastic inflation approach and we have studied 
the statistical distribution of curvature perturbations. They are very nearly Gaussian 
distributed: unlike other models where non scale-invariant fluctuations are produced, 
the feature in the spectrum of curvature perturbations does not imply a deviation 
from the Gaussian statistics. The class of scalar-tensor theory of gravity analysed is 
just one of the possible alternative theories to General Relativity that is consistent 
with observations. The results obtained show that, if gravity is given by one of this 
theories, the spectrum of curvature perturbations originated during chaotic inflation 
may be quite different from the usual one and may help to reconcile CDM models 

with large-scale observations. 
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FIG. 1. Power spectrum of C for the case (3 = 4 and 7 = 0.38. 

FIG. 2. Power spectrum of C for the case (i = 3 and 7 = 0.43. 


26 


ln(1/A(Mpc)) 



Fig. i 



((odiM)y/i)u| 



Fig. 2. 



